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ABSTRACT 

An iterative method for reconstructing mass distribution in spiral galaxies using a thin disk approximation 
is developed. As an example, the method is applied to galaxy NGC 4736; its rotation curve does not allow one 
to employ a model with a massive spherical halo. We find a global mass distribution in this galaxy (without 
non-baryonic dark matter) that agrees perfectly with the high resolution rotation curve of the galaxy. This mass 
distribution is consistent with the /-band luminosity profile with the mean mass-to-light ratio M/Li = 1.2, and 
also agrees with the amount of hydrogen observed in the outermost regions of the galaxy. We predict the total 
mass of the galaxy to be only 3.43 x 10 10 M Q . It is very close to the value predicted by the modified gravity 
models and much less than the currently accepted value of 5.0 x 10 10 M Q (with w 70% of the mass in a dark 
matter halo). 



Subject headings: 



Introduction 



Determination of dynamical masses of spiral galaxies is 
decisive step in establishing empirically the clustering 
^^cale of dark matter. The Cold Dark Matter model pre- 
F^licts a hierarchy of scales, from a dwarf sub-halo scale of 
of solar masses up to a galaxy cluster size. There is 
Qi strong evidence from the X-ray observations of galaxy 
^Slusters, and from detecting da rk matter with gravita- 
tional lensing Massev et al.ll2007 | for a 1 Mpc scale, which 
• is, the galaxy cluster scale. At the galactic size scale (of 
m £$he order of 10 — 50kpc), there are contradictory results. 
)ark matter was not obs erved in some elliptical galaxies 
2003l | . As for the spiral galaxies, it 



iRomanowsky et al 



as pointed out that the amount of dark matter in the 
disk of the Milky Way gal axy near the Sun is insignificant 
Kuiiken fc Gilmorelll989l ]. Moreover, mass estimates of 



spiral galaxies are model-dependent Sofue fc RubinlEoOl 

In the following we propose an Iterative Spectral 
Method to reconstruct mass distribution in the spiral 
galaxies. The method is bas ed on the thin disk model 
Toomrel[l96i lNordsiecld[l973] which offers the possibility 



to obtain mass density distribution fro m the observed rota- 



tion curve 



Binnev fc Tremaindll987l |. With the Iterative 
Spectral Method this reconstruction is unique, provided 
there is observational data about the hydrogen density 
in the regions outside of rotation velocity measurements. 
Without such data one can also extrapolate hydrogen 
measurements by assuming its falloff profile that would be 
physically acceptable. One should also take into account 
the abundance of helium which is proportional to that of 
hydrogen. In this way we obtain a complete set of data in 
the complementary regions of a given galaxy. Complete 



data are required to eliminate cutoff errors that normally 
result from cutting-off integration in an integral that re- 
lates rotational velocity to disk surface mass density at the 
radius where measurements of a rotation curve end. We 
shall call the radius the 'cutoff radius' and denote it by R. 

In fact, in early days of galactic research the thin 
disk model was applied to infer the galactic masses 
from the measured rotation curves. However, large er- 
rors of measurements of rotational velocities practically 
made the surface mass de nsity reconstruction impossible 
Binnev &: Tremainei 1987 1. To avoid uncertainties in the 



1 Jagellonian University, Institute of Physics, Reymonta 4, 30-059 
Krakow, Poland. 



thin disk mass reconstruction it became customary to a 
priori divide spiral galaxies into three major components: 
bulge, disk and spherical halo. Then a several parameter 
fit was made determining the properties of those compo- 
nents. Although this approach is physically appealing, 
the amount of arbitrariness in such a division is unaccept- 
able. In particular, the most massive component of this 
construction - the spherical halo - stands on the weakest 
empirical footing. 

In this paper we address the problem of determination 
of masses of spiral galaxies with rotation curves violat- 
ing the sphericity condition at large radii. One expects 
that the disk model should approximate mass distribution 
in such galaxies better than models assuming a massive 
spherical halo. The reconstructed surface mass density 
found for galaxy NGC 4736 approaches the hydrogen den- 
sity at the outskirts of the galaxy, the rotation curve is 
perfectly matched in the measured range, and its contin- 
uation is predicted further out into the Keplerian regime. 
Thus the disk model performs much better in this case 
than models including the spherical dark halo. 

To see why the assumption of the existence of spherical 
halo seems to be unjustified for a class of spiral galaxies, 
we consider as an example the galaxy NGC 4736. It was 



found (jKentl [1987[) that most of the mass (68%) in this 
galaxy is concentrated in a spherical dark halo, mainly 
at larger radii. At these distances surfaces of constant 
gravitational potential are almost spherical. Therefore one 
would expect in this region the rotation curve predicted by 
the Kent's analysis to satisfy the sphericity condition. For 
a spherically symmetric mass distribution the Keplerian 
mass function, Mk(p) = G^ 1 v 2 (p)p, corresponding to ro- 
tational velocity v(p), should be a nondecreasing function 
of the radial distance from the center /?, that is, 

v 2 (pi)pi < v 2 (p 2 )p2, if only p 1 < p 2 . 

The experimental rotation curve of NGC 4736 violates 
this condition for p > 6 kpc, see figure [TJ whereas the 
best-fitting solution for the rotation curve found by Kent 
satisfies this condition. At large radii the best-fitting 
rotational velocity is also greater from the observed by 
a factor of ps 1.2. The total mass should be therefore 
re 5.9 x 1O 1O M /(1.2) 2 = 4.1 x 1O 1 O M rather than 
5.9 x 1O 1O M as found in [Kentlll987j (at the currently 
accepted distance 5.1 Mpc to NGC 4736 the mass is even 
less, pa 3.5 x 10 10 M Q ). This shows that the mass distri- 
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approach was presented in [Kostovll2006] . where rotation 
curves were assumed Keplerian in the unobserved regions. 
This procedure is still unsatisfactory, as velocities of test 
particles on circular orbits encircling a disk-like object 
need not be Keplerian. Besides no comparison with the 
observed mass density has been done. 

2. Iterative Spectral Method 

Let vr = v(R) denote the last measured value of rota- 
tional velocity at cutoff radius R. We define also a dimen- 
sionless function u{x) = v(xR)/vr for describing rotation 
curve, where x = p/R measures distances in units of the 
cutoff radius (we work in cylindrical coordinates (/?, 4>, z)). 

Below we shall derive equations that are fundamental 
for the iteration method we present later on. 

2.1. Basic equations of disk model 

Let the z = plane represent the galactic disc. It is 
assumed that most of galactic mass is concentrated in its 
vicinity. A z-symmetric and cylindrically symmetric vac- 
uum gravitational potential is given by 

CO 

$(p,z) = -2ttv% J a(cj)J (lj^) e-^du. 
o 

A discontinuity of m ^ e direction normal to the 

z = plane is interpreted as a surface mass density a(p) 
of the galactic disk 

oo 

°{P) = || / «>&(u)J dw, (2-1) 



thus <?{lo) is the Fourier-Bessel transform of a(p) 

oo 

GR f 

ct(uj) = — 2~ / x a(Rx)Jo(u> x)dx. (2-2) 
v r J 



Fig. 1. — Negative sphericity test for spiral galaxy NGC 
4736. The figure shows experimental rotation curve (solid 
line) and the corresponding Keplerian mass function de- 
fined by M K (p) = G- 1 pv 2 (p) (dashed line). The Ke- 
plerian mass function is not increasing at all distances, 
therefore a massive spherical halo in the approximation of 
circular orbits is excluded and should be replaced with an 
extended flattened component. 

bution in the galaxy cannot be dominated by a spherical 
component and that it would be more appropriate to as- 
sume more extended disk components in similar galaxies 
in place of massive spherical halos. 



Modelling of galactic rotation curves usually assumes 
that orbits of stars are circular and concentric, then 
y2 (p)\z=o = /° 5 p $ (/°'°)' therefore, 



u (x) — 27r x J lu (t{oj) Ji(ujx)duj. 



By substituting (|2-2|) to (|2-3| one can show thalQ 

. . n ! / W dx 

p - X 



(2-3) 



v\p)=AGp V\ a( X ) 



L We used the following identities Ryzhik & Gradshtein 1951]: 



Th e global disk model has been used recently [Takamiva fc 
2000]. however the procedure of cutting-off integrals ap- 
plied therein lead to overestimation of surface density at 
larger radii for galaxies with declining rotation curves. For 
example, surface density in the galaxy NGC 4736 could 
be reliably determined only within a radius of less than 
6 kpc. But the experimental data available for this galaxy 
enable one to find a much better determination of mass 
distribution as far as 16 kpc from the center. An improved 
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7T 1 [K(x/0-E(x/0]x-\ 0<x<£, 



where the complete elliptic functions are defined by 

it/2 it/2 

*(*)=/ A I- ■ « 

J vl-K sin 
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5(k) = J tyy/ 




1 — k, sin 



<r{x) 



X 2 E 



(i) 



p(x 2 - P 2 ) 



(2-4) 



where "P stands for the "principal value integral". The 
simple pole x — P dominates the logarithmic divergence in 
K(p/x), thus the principal value integral is easily tractable 
numerically and can be calculated with high accuracy. The 
inverse of the above integral read^l 



v 2 ( X ) 



P X 




E 



(i) 



*G V \] V {X) X 2 -P 2 



dx 



E 



(i) 



X P 2 - X 2 



dx 



(2-5) 



Unfortunately, this nice formula is not very useful for de- 
termining <j(p) directly from v(p), as the observed rotation 
curves have naturally a cut-off, and therefore the integra- 
tion cannot be carried outd For example, the closer p is 
to cutoff radius R, the more <j(p) is underestimated by 
cutting off the integration at R. In order to overcome this 
unpleasant feature of disk geometry one may use another, 
complementary set of the experimental data concerning 
mass distribution at p > R. A global mass density and 
the corresponding global rotation curve can be found it- 
eratively that would both satisfy the relations (|2-4p and 
(12-51) and be consistent with the observations. 



2.2. A series approximation of the surface density 

Below we shall find a first approximation to cr{p) di- 
rectly from the observed part of rotation curve. 

A function u 2 {x)/x defined on the unit interval x € 
(0,1), can be expressed in terms of a series of indepen- 
dent functions. In particular, if u 2 (x)/y/x is integrable 
over the interval, a complete set of orthogonal functions 
on x € (0, 1) can b e constructed from the Bessel function 
Ji |LebiedievlH95H . Then 



u 2 (x) 



k 



< x < 1, Jo(^fc) = 0, uo k > 0. 



(2-6) 



The summation is taken over all consecutive zeros u>k > 
of the Bessel function Jo. By multiplying both sides by 
x Ji(uj m x) and integrating over x € (0, 1), we find that 



2 f 2 Ji{uj k x) 
Pk, Pk = — u {X)— —. — -da;. (2-7) 

J J\\Wk) 



The amount of information encoded in the enumerable 
sequence is insufficient for determining the transform 
<7(lu) which is a continuum of numbers. Therefore, a(p) 



2 Equation 112-51 1 is deliberately presented in the form containing only 
squares of rotati onal velocity. Much more familiar is its equivalent 
form given in the lBinnev fc Tremaine] [1987H handbook, and therefore 
commonly used, however it has the deficiency of containing deriva- 
tives of v 2 (p) which are determined experimentally with very low 
accuracy. 

3 however, the relation serves as a consistency check of our final results 
we obtain iteratively later on 



cannot be determined accurately from the observed rota- 
tion curve alone, even in the region p < R. However, the 
elements of the sequence can be still linearly combined to 
approximate a(ui). Note, that every such approximation 
is equivalent to choosing a particular extrapolation of a 
rotation curve beyond its cutoff radius. 

There are some mathematical subtleties regarding this 
extrapolation. Let us assume first that we put v(p) — for 
p > R and v(p) — vr/2 at a single point p = R. Then we 
can use the representation f|2-6[) for x < 1, and for x > 1 
we set u — 0. From the inverse transform of (|2-3p . which 
reads 



<j(uj) = — I u 2 (x)Ji(uix)dx, 
2n 



we then obtain a continuous spectrurrjfl 



2tt ^ 

k 



OjJ (u)Ji{uJk) 



(2-8) 



(2-9) 



However, for handling the observational data more tractable 
are discrete spectra. For instance, by extending formally 
(|2-6j) to the region x > 1, we get 



k 



(2-10) 



Now, the corresponding u 2 (x), calculated from equation 
(|2-3[) , can behave quite odd for x > 1 and become even 
negative! Nevertheless, both (|2-9p and (|2-10[) . when sub- 
stituted to (|2-3|) , lead to the same u 2 (x) for any x £ (0, 1). 
Although the corresponding surface mass densities differ 
from each other also for p < R, the difference is, in prac- 
tice, negligible close to the centre. We may therefore ac- 
cept (|2-10[) as a first approximation of d{u>). By substitut- 
ing (pTTOl) to J2^T|) we obtain 



ip) 



k 



(2-11) 



Despite that close to R equation (|2-11|) is not a trustworthy 
approximation (the sum may become negative), it almost 
perfectly agrees (usually for p < 0.6R) with the surface 
densities corresponding to several model rotation curves 
that were cutting-off at different R. 

2.3. The iteration scheme 

From the above, it is evident that the total mass of the 
galactic disk within the cutoff radius cannot be determined 
from its rotation law only. However, if additional bounds 
for the surface mass density are known outside the cutoff 
radius, one may find a reliable global mass distribution by 
iterations. 

As before we denote by R the cutoff radius, to which 
a rotation curve has been measured. We denote this part 
of rotation curve by i>o- The distribution of matter, which 
extends far beyond the cutoff radius and which is mostly 
composed of neutral and molecular hydrogen, we denote 
by <Jh ■ If such data are not available, one may try to devise 
a physically plausible profile of baryonic matter outside R. 



4 Since Jo(w) ~ Jl(Vk){uk ~ w ) when u> ?a the components under 
the summation sign are continuous functions of u> 



At the very beginning we define u(x) = vq(Rx)/vo(R). 
We calculate from (|2-lll) an initial surface density Co where 
the corresponding spectral coefficients are calculated 
from (|2-7p (we put <tq and in place of a and &k , respec- 
tively). Let i?i be the first value of p where cto(p) = o'nip)- 
The first approximation a\ to the sought global surface 
density for p < R\ is defined by setting a± = an, and for 
p > R\ we set <j\ = cro. With this o\ we begin the first it- 
eration step (i = 1) which is described below for arbitrary 



The 



i iteration step By setting a — oi in equation 
(|2-4[) we calculate the corresponding approximation 
v — Vi to a global rotational velocity. Afterwards we 
find the resulting discrepancy function 8vf — v q — v\ 
defined for p < R. The function measures an ap- 
proximation error in the i th iteration step and is 
used also to determine a correction Sat to the ac- 
tual global surface mass density Oi. By substitut- 
ing 5vf (x R) I Svf (R) in place of u 2 {x) in equation 
(|2-7[) we find spectral coefficients 5a^ where now &k 

stands for Sa^ . Next we calculate for p < R the 
correction 5ai from equation 



8<Ji(p) 



5vf{R) ^ (i) 
fe 



P 

Uk R 



which is analogous to equation (|2-11[) . Let Ri+\ be 
the first value of p where (Ji{p)+5<Ji{p) = <jh{p)- The 
corrected global surface mass density c^+i is found 
by setting o-j+i = <n + 8a i for p < Ri+\ and o-j+i = 
a H for p > Ri+i. 

Now the i + 1 iteration step follows. We repeat this proce- 
dure until max(Sv 2 l ) in the n th iteration becomes accept- 
ably small. 

3. Application of the iteration method to the spi- 
ral galaxy NGC 4736 

As an example we choose the spiral galaxy NGC 4736. 
Its rotation curve excludes the massive spherical halo, 
therefore the global disk model should better approximate 
the gravitational field in this galaxy than models assuming 
such a halo. 

There is also the additional information about distri- 
bution of hydrogen observed in the remote region of this 
galaxy, where the rotation curve has not been determined. 
This provides the lower bound for the surface density. 
No other luminous matter has been observed there. The 
amount of hydrogen in this region is comparably small. 
Its presence allows us to reconstruct the global mass den- 
sity and predict the rotation curve beyond the presently 
measured range. 

3.1. Observational data 

A high res olution rotation c urve of the galaxy was 
adopted from [Sofue et al . 1999]. It was measured out 
to 7' which corresponds to cutoff radius R = 10.38 kpc 
(the distance of the galax y is D = 5.1 Mpc and inclina- 



tion angle i = 35° as in [Sofue et al.l fl999]). The rota- 



Sv — ±llkm/s (attained at p — 5') and the average ve- 
locity is 168 km/s. Thus the total dynamical mass should 
be determinable with the accuracy better than w 10%. 
We use also the 50" resolution neutral hydroge n surface 



We use also tne ou resolution neutral nydroge n surlaci 
density published in [Mulder fe van Driell 1993| (inclina- 



tion was corrected from 40° to 35°), which is consistent 



with [Bosma et al.lll977j |. A /3-band surface brightness A/3, 
where (3 = I, V or B is calculated from 



K 



L 



0,& 



pc^ 



cos (i) 



tan (arcsec) 



_ 10 0.4(M 3 , o -^-5)_ 



It is done using the respective luminosity profiles, cor- 
rected for galactic extinc tions nn L™„ a „ s „2 l < taken from 



i ' ML ckuc j ' 

Munoz-Tunon et al.lll989l | . The refer ence absolute magni 
tudes for the Sun a re: Mj, P) = 4.19 Ivan den Bosch at al 
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4.83 



5.43 |Carignanl[l985 



Takamiva fe SofudbOOOj and M 



'B,0 



3.2. 



Determining the global surface density and 
the global rotation curve 



Below we shall apply the iteration scheme described in 
section 12.31 We use the same notational conventions as 
introduced there. 

A third order interpolating function u 2 (x) = Vq(Rx) jv 2 R 
(vr = v (R) — 125.6 km/s) is used to find the spec- 
tral coefficients The pu coefficients in (|2-7[) fall 
off quickly to zero, eg. Pi/pi2b ~ 3 x 10~ 6 , and the 
functional sequence of partial sums in (|2-11|) converges 
quickly. The summation can thus be terminated at 
some finite k (we put k = 144). This truncated sum 
is taken as the surface mass density (Jo(p). Next, we 
set ai(p) = (T (p) for p < Ri, <7i(p) = <r H {p) for 
Ri < p < Rh and a\{p) = elsewhere, where (Jh{p) 
is the azimuthally averaged neutral hydrogen surface den- 
sity, Ri was defined in sec. 12.31 and Rh = 16.06 kpc > R, 
cf . line '(7i ' in figure [2] The integrated mass of o\ (p) is 



NGC 4736 




tional velocity has been measured with the maximum error 



p [kpc] 



Fig. 2. — A fragment of surface mass density of galaxy 
NGC 4736 obtained in consecutive iterations (solid lines 
<7i, 02 j 03 5 04) and the observed column mass density of 
Hi (erff - solid circles). 

Mi « 3.144 x 10 10 M Q and should not exceed the total 



mass Mtot more than < Mtot ~ Mi < 8M\, where 
SM 1 < 2Rv (R)Me,x\v - v x \/G w 5 x 10 9 M Q w O.I6M1. 
The corresponding first approximation «i (p) to the global 
rotation curve agrees with vo(p) very well for the small 
radii, where &i(p) is reliable. However, there is still a 
discrepancy at greater radii, cf. line v\ in figure [3l In 
order to reduce it, we find a correction Sai(p) to cri(p) in 
the region < p < R. The needed spectral coefficients 
cffc were determined by finding a least square fit to a 
first order interpolation of Svf(p), with xJi(cJkx) used 
as basis functions!! Having found Sai(p) we calculate 
the corrected surface density in the second approxima- 
tion 02 (p) = fi (p) + <5fi (p) ! with the reservation we set 
02 (p) = ch{p) for P> Ri and o"2(p) = for p > (the 
line '172' in figure [5]). Now, the corresponding global ro- 
tation curve i>2(p) agrees with v (p) almost perfectly (the 
line 'i>2 ' in figure [3]) , but it still needs to be corrected a 
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Fig. 3. — vi,V2,V3,V4 are rotational velocities calculated 
from equation ()2-4|l for global surface densities obtained in 
the four consecutive steps of iteration. The thin line V4 was 
shown also outside the observational range (experimental 
points for NGC 4736 are indicated by circles). 

little in the vicinity of, R, since Max(l — ^2/^0) ~ 0.02. 
The corresponding mass is M2 = 3.385 x 10 10 M© and is 
underestimated, but Mtot — M2 < 8M2, where 8M2 < 
2Rv Q (R)Ma,x\v Q - v 2 \/G « 1 x 1O 9 M » 0.03M 2 . Simi- 
larly, we carry out the third step, to find 03 (p) and v^{p). 
The integrated mass in this step is M 3 = 3.419 x 1O 1O M 
and SM 3 w 3.3 x 10 8 M Q w O.OIM3. The error has not 
changed substantially compared with 8M2, therefore we 
terminate iterations after the fourth step. In this way we 
obtain the global surface density 04 (p) shown in figure 
IU It gives the total mass M TO t = 3.43 x 10 10 Af Q (it 
includes the observed hydrogen Hi out to Rh), and the 
corresponding rotation curve v±{p) - its predicted contin- 
uation beyond the cutoff radius can be tested with future 
observations (see figure [5]). 



5 It is known from the theory of generalized Fourier series, that in the 
theoretical limit of infinite number of terms in the sum H2-61 . the 
values of a k determined by the least square method would be equal 
to those calculated from equation 112-71 1 



Fig. 4. — Global surface mass density for NGC 4736 found 
in the fourth iteration (thick solid line 04). This mass dis- 
tribution leads to rotation velocity that agrees perfectly 
with the observed rotation c urve. For comparison, t he 
surface mass density found in Takamiva k, Sofue1l2000| | is 
also shown (the thin line s). The observed density pro- 
file of hydrogen Hi is marked with solid circles. Luminos- 
ity profiles A(p) (at inclination angle 35°) in the i3-band 
(dotted triangles), the I/-band (dotted squares) and the 
/-ba nd (dotted circles) were cor rected for galactic extinc- 



tion |Munoz-Tunon et all [1989} . The /-band luminosity 



was corrected also for the internal extinction (marked with 
stars). 

3.3. Discussion 

In the previous section we have obtained the surface 
density in galaxy NGC 4736 that is consistent with dis- 
tribution of Hi beyond the cutoff radius R, and perfectly 
agrees with the observed rotation curve. Since Hi mass 
outside R is only 0.003Mtot, the results (e.g density pro- 
file) would not be affected significantly if there was an 
amount of baryonic matter outside R comparable to that 
of Hi. The global surface mass density of galaxy NGC 4736 
found in the fourth iteration is shown in figure [4] (the '(74' 
line). The corresponding global rotation curve calculated 
from equation (|2-4p agrees perfectly with the experimental 
rotation curve vq, see figure [5]) • 

The V and B-band luminosities are smaller than the 
/ band luminosity Li, cf. table Q] (Ly ~ 0.7Lj and 
Lb ~ 0.5-Lj), which may indicate the presence of dust. 



Mtot 


3.43 x 10 


iO M 


M m 


6.00 x 10 8 M Q = 


>■ 0.02A/ TOT 


Li 


2.87 x 10 10 L 7iQ = 


> M/Lj = 1.2 


L v 


2.14 x \Q w L v .q = 


> M/L v = 1-6 


L B 


1.54 x 1O 1o Lb' = 


> Mj L B = 2.2 



Table 1: Physical parameters obtained with the help of the 
Iterative Spectral Method for spiral galaxy NGC 4736. (at 
the distance 5.1[Mpc] and inclination angle 35°). 

The mass-to-light ratio profile (surface density to surface 
luminosity ratio) is shown in figure [6] From the three 
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NGC 4736 



1.35' 2.7' 4.04' 5.39' 6.74' 8.09' 9.44' 10.79' 




small: M/Lj w 1.2. Also its distribution, i.e. den- 
sity to luminosity ratio which is of the order of 1 - 2 
(cf. figure O is typical to stars of slightly sub-solar mass. 
The total galaxy mass is surprisingly small. It is com- 
parable to the mass predicted in the framework of the 
MOND model, which gives 3.21 ± 0.09 x 10 10 A/ Q , or in 
the framework of the metric skew tensor gravity, which 



Brownstein fc Moff at 2006]. 



gives 3.15 ± 0.08 x 1O 1U M 
We obtain M T ot = 3.43 x 10 10 A/ o for the same rotation 
curve in the framework of Newtonian physics. 

In figure [7] is shown the mass function 

M 4 (p)/M TOT , M 4 (p) = 2ir f <r 4 (£)£d£. 

Jo 

It is shown together with the Keplerian mass function 
Vq(p)p/(GMtot), corresponding to the observed rotation 
curve vo(p). The 'odd' behaviour of the latter function was 



Fig. 5. — The observed rotation curve of galaxy NGC 
4736 (circles) and the rotation curve (solid thick line v 4 ) 
corresponding to surface mass density <t 4 found in the 
fourth iteration and calculated from equation (|2-4|) . The 
solid thin line Vk is the corresponding Keplerian asymptote 
^GMtot/p where Mtot ~ 3.43 x lO lo A/ is the total 
galaxy mass. The solid thin line v m is the rotational veloc- 
ity which one would obtain from a spherically symmetric 
mass distribution with the same mass function as that for 
NGC 4736. 
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Fig. 6. — The predicted ratios of mass-luminosity sur- 
face densities for galaxy NGC 4736 with luminosities cor- 
rected for galactic extinction in B-band (triangles), V- 
band (squares), /-band (circles) and corrected additionally 
for internal extinction /-band (marked with stars). 



bands, the /-band should be regarded as the most reli- 
able tracer of the luminous matter, as it most accurately 
maps the amount of the radiated energy. Luminosity in 
V- and _B-bands is reduced by the presence of dust, it 
also does not include the radiation from low mass stars. 
For galaxy NGC 4736 the /-band mass-to-light ratio is 




Fig. 7. — Mass function for NGC 4736 normalized with 
respect to total mass Mtot = 3.43 x 10 10 M Q (thick solid 
line) compared with Keplerian mass function m(p) (my 
- the thin solid line) defined as pvq(p)/(GMtot), where 
vq(p) is the observed rotation curve. 

our argument against applying dark halo dominated mod- 
els to galaxy NGC 4736. Such a model used in [KentHl987j 
predicted 5.0 x 1O 1O M for the total galaxy mass at a dis- 
tance 5.1 Mpc. The best fitting curve in that model gives 
rotational velocity at the cutoff radius 1.2 times greater 
than the observed value. As this leads to overestimation 
of the total mass by a factor of some (1.2) 2 , one should ex- 
pect rather 5.0 x 10 10 A/ Q /(1.2) 2 « 3.47 x 1O 1O M for the 
total mass, which is very close to our result obtained in the 
global disk approximation. We note also a related and triv- 
ial observation: Keplerian rotation law ^/ M/^p)G j p for a 
spherically symmetric matter distribution with the mass 
function numerically equal to Af 4 (p), gives smaller rota- 
tion velocities at larger radii than the disk model. Also 
more matter would be required by a spherical model to 
account for observations, see figure [5l 



4. Conclusions 

The comparison of mass functions and rotation laws at 
the end of the previous section, illustrates the fact that 
the models with flattened mass distributions are more effi- 
cient than the commonly used models assuming spherical 
halo. The former are better in accounting both for high 
rotational velocities as well as for low scale structure of 
rotation curves and with noticeably less amount of mat- 
ter than the latter (the relation between rotation and mass 
distribution in the disk model is very sensitive for gradients 
of a rotation curve). The use of the disk model is justified 
for galaxies with rotation curves violating the sphericity 
condition. This is necessary (although not sufficient) con- 
dition for a spherical mass distribution^ 

Rotation of the spiral galaxy NGC 4736 can be fully un- 
derstood in the framework of Newtonian physics. We have 
found a mass distribution in the galaxy that agrees per- 
fectly with its high-resolution rotation curve, agrees with 
the /-band luminosity distribution giving low mass-to- light 
ratio of 1.2 in this band at total mass of 3.43 x 1O 1O M0, 
and is consistent with the amount of Hi observed in the 
remote parts of the galaxy, leaving not much room (if any) 
for dark matter. Remarkably, we have achieved this con- 
sistency without invoking the hypothesis of a massive dark 
halo nor using modified gravities. 

There exist a class of spiral galaxies, similar to NGC 
4736, that are not dominated by spherical mass distribu- 
tion at larger radii. Most importantly, in this region ro- 
tation curves should be reconstructed accurately in order 
not to overestimate the mass distribution. For a given ro- 
tation curve it can be easily determined whether or not a 
spherical halo may be allowed at large radii by examining 
the Keplerian mass function corresponding to the rotation 
curve (the so called sphericity test). 

By using complementary information of mass distribu- 
tion, independent of rotation curve, we overcame the cutoff 
problem for the disk model, that for a given rotation curve, 
a mass distribution could not be found uniquely as it was 
dependent on the arbitrary extrapolation of the rotation 
curve. 



Kent, S. M. 1987, AJ, 93, 816 

Kostov, V. 2006, preprint ( |astro-ph/0604395fr d) 

Kuijken, K., and Gilmore, G., 1989 MNRAS, 239, 605, 651 

Lcbiediev N. N., 1953, Special functions and their appli- 
cations (Moscow) 

Massey, R. 2007, Nature, 445, 286 

Mulder, P. S., and van Driel, W., 1993 A&A, 272, 63 

Munoz-Tunon, C, Prieto, M., Beckman, J., and Cepa, J., 
1989, Ap&SS156, 3001 

Nordsieck, K. H. 1973, ApJ, 184, 735 

Romanowsky, A. J. 2003, Science, 301, 1696 

Ryzhik, I. M., and Gradshtein, I. S. 1951, Tables of Inte- 
grals, sums, series and products (Moscow, Leningrad) 

Sofue, Y., and Rubin, V. 2001 ARA&A39, 137 

Sofuc, Y., Tutui, Y., Honma, M., Tomita, A., Takamiya, 
T., Koda, J., and Takeda, Y. 1999, ApJ, 523, 136. 

Takamiya, T., and Sofue, Y. 2000, ApJ, 534, 670 

Toomre, A. 1963, ApJ, 138, 385 



One of the authors (MK) was supported by the Polish 
Ministry of Science and Higher Education, grant 1 P03D 
005 28. 



REFERENCES 

Binney, J., and Trcmainc, S. 1987, Galactic dynamics 
(Princeton Univ. Press) 

Dutton, A. A., van den Bosch, F. C, Dekel, A., and 
Courteau, S. 2006, preprint ( |astro-ph/0'604 553vl) 

Bosma, A., van der Hulst, J. M., and Sullivan III, W. T. 
1977, A&A, 57, 373 

Brownstein, J. R., and Moffat, J. W. 2006, ApJ, 636, 721 

Carignan, C. 1985, ApJS, 58, 107 



6 A thin disk with surface density (t(x) = (1 + x 2 ) 3 / 2 / (2tt) rotates 
with velocity v(x) = x(l + x 2 )~ 3//4 . This rotation low satisfies the m , T 

sphericity condition: (xv 2 (x))' = 3x 2 (l + x 2 )"5/ 2 > 0. Thls 2 - column P re P" nt was prepared with the AAS IATeK macros 

v5.2. 



